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Starting from hyperbolic dispersion relations for the invariant amplitudes of pion-nucleon scat- 
tering together with crossing symmetry and unitarity, one can derive a closed system of integral 
equations for the partial waves of both the ^-channel (jiN — » 7lN) and the t -channel (%% — » NN) 
reaction, called Roy-Steiner equations. After giving a brief overview of the Roy-Steiner system 
for tcN scattering, we demonstrate that the solution of the ^-channel subsystem, which represents 
the first step in solving the full system, can be achieved by means of Muskhelishvili-Omnes tech- 
niques. In particular, we present results for the P- waves featuring in the dispersive analysis of the 
electromagnetic form factors of the nucleon. 
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1. Introducing Roy-Steiner equations for 7iN scattering 



Partial-wave dispersion relations (PWDRs) together with unitarity and crossing symmetry as 
well as isospin and chiral symmetry (i.e. all available symmetry constraints) have repeatedly proven 
to be a powerful tool for studying processes at low energies with high precision For 7iN scat- 

tering the (unsubtracted) hyperbolic dispersion relations (HDRs) for the usual Lorentz-invariant 
amplitudes read [Q] (using the notation of [^], see [f7|] for more details) 
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and similarly for A - , B~, and N~, where A/^ are the nucleon pole terms and the "external" (un- 
primed) and "internal" (primed) kinematics are related by real hyperbola parameters a and b (as 
well as via s + 1 + u = 2(m 2 + M 2 ) = s r + t r + u r ), so that HDRs allow for the combination of all 
physical regions, which is known to be crucial for a reliable continuation into the subthreshold 
region and hence for an accurate determination of the TtN <7-term. Furthermore, the imaginary 
parts are only needed in regions where the corresponding partial-wave decompositions converge 
and the range of convergence can be maximized by tuning the free hyperbola parameter a. While 
the ^-channel integrals start at the threshold s + = W 2 = (m + M K ) 2 , the ^-channel contributes al- 
ready above the pseudothreshold t % = AM 2 far below the threshold f# = 4m 2 . Depending on the 
asymptotic behavior of the imaginary parts, in principle it could be necessary to subtract the HDRs 
to ensure the convergence of the integrals, thereby parameterizing high-energy information with 
polynomials containing a priori unknown subtraction constants. However, (additional) subtrac- 
tions may also be introduced to lessen the dependence of the low-energy solution on high-energy 
input; the corresponding subtraction parameters then obey respective sum rules. For kN scatter- 
ing it proves particularly useful to subtract at the subthreshold point {s = u,t = 0), as this pre- 
serves the s <-> u crossing symmetry (which can be made explicit in terms of the crossing variable 
V = (s- u)/(4m) via D ± (v,t) = A ± + vB ± = ±D ± (-v,t)). This is especially favorable for the 
^-channel subproblem and facilitates matching to chiral perturbation theory ^] to determine the 
subtraction constants, which thus can be identified with the subthreshold expansion parameters. 1 
In addition to the presentation in [^], we also introduce a (partial) third subtraction, which is related 
to the parameters and a^ x of the subthreshold expansions (with d^ n = d^ n for all n > 0) 

A + (v,0 = — + J++J+^ + a+v 2 + ^(v 4 ,v 2 /^ 2 ) , A~(v,f) = « 00 v + a " 1 vr + ^(v 3 ,v/ 2 ) . 

(1.2) 



^or the PWDRs of kk scattering, called Roy equations [|iq|], an analogous matching procedure for the kk scattering 
lengths as pertinent subtraction parameters has been conducted in fH. In contrast to Tin scattering, the 7iN scattering 
lengths can be extracted with high accuracy from hadronic-atom data [|l2| |l3|] and may thus serve as additional con- 
straints on the subtraction constants in the Roy-Steiner system. 
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In order to derive the partial-wave HDRs, called Roy-Steiner (RS) equations, one needs to 
expand the s- and /-channel imaginary parts in dlTT] ) into the respective partial waves and subse- 
quently project the full expanded equations onto either s- or /-channel partial waves; the resulting 
sets of integral equations together with the respective partial-wave unitarity relations then form the 
s- and /-channel RS subsystems. According to [||], the (unsubtracted) ^-channel RS equations read 
(based on the MacDowell symmetry fl +l) _{W) = -f l+ {-W) for alH > @) 

oo 

f}+ (W) = N' e+ (W) + X - JdtX {Gu(W/) Jmf J + (t 1 ) + H U (W, t') lmf J _ (t') } 
u 3 

oo 

+ 1 l m ' £ {Kl ff {W,W')lmf e+ {W t )+ , (1.3) 

w + 

where due to G-parity only even/odd / contribute for isospin / = +/—, respectively, and the partial- 
wave projections of the pole terms as well as the (lowest) kernels are analytically known, the latter 
including in particular the Cauchy kernel: K T U , (W, W') = 8w / (W f — W) + ... . The ^-channel / = ± 
partial waves are intertwined by the usual unitarity relations, which are diagonal in the ^-channel 
isospin basis I s E {1/2, 3/2} only. Once the /-channel partial waves are known, the structure of the 
^-channel RS subsystem is therefore similar to the Tin Roy equations, cf. [[[[]. As shown in [0], the 
corresponding (unsubtracted) /-channel RS equations are given by 

oo 

fi(t)=Ni(t) + ±- /dW'£{G//(^W0lm^ + (W0 + Gy/^-W0lmyf m) _(W0} 

oo 

+ ^/^T{^'( f ' ?, )Im4'(^) + ^(^0lm/'(0} (1-4) 

and similarly for the f_ except for the fact that these do not receive contributions from the f+. 
Here, only even or odd /' couple to even or odd / (corresponding to /-channel isospin I t = or 
I t = 1), respectively, and Kj Jf (as well as the analogous Kj Jf for the f J _) contains the Cauchy kernel. 
Moreover, it turns out that only higher /-channel partial waves contribute to lower ones. Assuming 
Mandelstam analyticity, the equations (|L4|) are valid for ^ft <E [2M n , 2.00 GeV] using a = -2.7 1M^, 



whereas ([TJ) holds for W e [m + Af^, 1.38 GeV] using a = -23.19M^. The /-channel unitarity 



relations are diagonal in I t and only linear in the f± (below the first inelastic threshold /i ne i) 
Im/iW = o?$(t)yfi(t)0(t-t x ) , o?t){t) = sm8j'(t)e iS j'^ , o?(t) = ^T^A , 

from which one can infer Watson's final state interaction theorem Jl5[ ] stating that (in the "elastic" 
region) the phase of f± is given by the phase Sj of the respective kk scattering partial wave tj . 

Due to the simpler recoupling scheme for the f±, the /-channel RS subsystem can be recast as 
a (single-channel) Muskhelishvili-Omnes (MO) problem jl6| ] with a finite matching point t m 



for and the linear combinations T\t) = my/j /{J +\)fi (/) - f{ (/) with T J {t N ) = for all 

/ > 1 of the generic form (the details are given in ^) 
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Figure 1: Flowchart of the solution strategy for the Roy-Steiner system for nN scattering. The third step 
consist in the self-consistent iteration (denoted by thick arrows) of the preceding steps until convergence. 



where the inhomogeneities A(t) subsume the nucleon pole terms, all ^-channel integrals, and the 
higher /-channel partial waves. For t % <t < t m , solving for \f(t)\ only according to Watson's 
theorem requires 8(t) for t n <t <t m and \mf(t) for / > t m . Introducing n > 1 subtractions does 
not change the general structure of the RS/MO system, e.g. the P-waves are given by 

00 00 
r W -A rW + - J tln . l{t ,_ tN){tl _ t y /-W-A-W +~J t , n{t ,_ t) , 

demonstrating that F J and hence f{ is effectively subtracted by one power less than f_, which 
motivates the additional (partial) third subtraction in A ± , cf. (|L2|), that affects solely the f+. 

The solution strategy for the full RS system in the low-energy (or even subthreshold/pseudo- 
physical) regions, where only the lowest partial waves are relevant and inelastic contributions may 
be (approximately) neglected, is shown in Fig. [j]; see ^ for more details. 

2. The /-channel Muskhelishvili-Omnes problem: P-wave solutions 

As the first step in the numerical solution of the full RS system, we check the consistency of 



our /-channel MO solutions with the results of the KH80 analysis [|17|], which are still used nowa- 
days although no thorough error estimates are given (and despite the availability of more modern 
experimental data). Here, we present results for the P- waves in the (elastic) single-channel approx- 
imation of the MO problem, which is well justified for the P- and higher partial waves, whereas the 
S- wave requires a two-channel description including KK intermediate states as described in [^]. To 
produce the results (that will also serve as input for the solution of the ^-channel RS subsystem, cf. 
Fig. [j]) partly shown in Fig. 0, we have used as input nn phase shifts from [|18|], ^-channel partial 
waves (/ < 4) from SAID JEJ| for W < 2.5 GeV, and above the Regge model of [gOj]. To facilitate 
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Figure 2: ^-subtracted MO solutions for the P-wave moduli. 



— 1-sub 

2- sub(a — > oo 

— 2-sub __ 

-- 2-sub(VU = 1.1 GeV) 

3- sub(a — > oo 
— 3-sub 

3-sub(V^ = 1.1 GeV) 
x KH80 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Vt I GeV 



Figure 3: Two-pion-continuum contribution to lmG v E and ImG]^ 



comparison with the results of KH80, we use the respective subthreshold parameter values and a 
7iN coupling of g 2 / (47r) = 14.28 [TtJ] (as starting point, the final values will result from the 



iteration procedure, cf. Fig. [j]). 2 Moreover, KH80 uses different types of dispersion relations, in 
particular so-called fixed-/ ones, which can be emulated (up to the /-channel contributions that are 
not present at all in the fixed-/ case) by taking the "fixed-/ limit" \a\ —> oo. As argued in [0], all 
/-channel input above ^Jt^ = 0.98 GeV is set to zero, which forces the MO solutions to match zero 
at / = / m . While Fig. |2| displays the results for \a\ — » oo, investigating the effect of using a different 
(i.e. higher) matching point leads to the same conclusion: with increasing number of subtractions, 
thus lowering the dependence on the high-energy input by introducing more subthreshold param- 
eter contributions as subtraction polynomials, the solutions show a nice convergence pattern both 
in general (proving the internal consistency and numerical stability of our RS/MO framework) and 
in particular towards the KH80 results (being consistent with relying on KH80 values for g and 
the subtraction parameters). The P- waves feature prominently in the dispersive analysis of the 
nucleon electromagnetic form factors, see e.g. and references therein, and in Fig. ^ we illus- 
trate the effects on the spectral functions (by approximating the vector pion form factor F% via a 



2 Modern analyses yield significant smaller values for the izN coupling, cf. e.g. g 2 / (4k) = 13.7 ±0.2 of [O. 
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twice- subtracted Omnes representation, cf. [[J) 

toG£(0 = ^(^(0)71(00^-^) , hnGUt) = t -^(F^t)YfUt)e(t-t 7: ) . 

We are confident that a self-consistent iteration procedure between the solutions for the s- and 
/-channel eventually will yield a consistent and precise description (including error estimates) of 
the low-energy 7iN scattering amplitude in all kinematical channels. 
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